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Abstract 

Calculations are performed of the cooling of neutron stars with standard and enhanced 
neutrino energy losses in the presence of neutron and proton superfluidities in the stellar 
cores. The effects of superfluidity on the heat capacity and the neutrino luminosity produced 
by the direct and modified Urea processes and by the neutrino bremsstrahlung emission in 
nucleon-nucleon collisions are taken into account. The constraints are analyzed on the critical 
temperatures T cn and T cp of the transition of neutrons and protons to the superfluid state, 
which can be obtained from a comparison of observational data on thermal radiation from 
neutron stars with theoretical cooling curves for the standard and enhanced neutrino energy 
losses. Possible ranges of T cn and T cp for the pulsar Geminga are discussed. 

1 INTRODUCTION 



Neutron stars (NSs) are unique astrophysical objects. The density in the NS central layers is 
several times higher than the standard nuclear density po = 2.8 x 10 14 g-cm -3 . Such a dense 
matter cannot be analyzed in laboratories, but it can be studied by astrophysical methods, in 
particular, by comparing theoretical calculations of the cooling of NSs with observational data on 
their thermal emission. 

The nuclear composition of the NS central layers is not known with certainty (Shapiro and 
Teukolsky 1983). We will restrict our consideration to the simplest model in which matter is 
composed of neutrons (n), protons (p), and electrons (e). The cooling rate of a not-too-old star 
(with an age t < 10 5 yr) is determined by the neutrino luminosity of its inner layers. One can 
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distinguish two types of NSs, with regard to the efficiency of the neutrino energy losses. Stars 
with a relatively low (standard) energy loss rate (see, e.g., Friman and Maxwell 1979; Yakovlev 
and Levenfish 1995) caused by the modified Urea processes 

n + N~^p + N + e + u e , p + e + N^n + N + u e , (1) 

and by the bremsstrahlung radiation of neutrino pairs in nucleon-nucleon (N) collisions 

N + N-*N + N + is + v. (2) 

constitute the first type of NS. 

A buffer nucleon iV in reaction ([D) is needed to satisfy momentum conservation in dense matter 
with relatively low number densities of e and p. To avoid misunderstandings, note that we are 
dealing here with highly degenerate neutrons and protons. The neutron branch of the reaction 
(N = n) is well known; the proton branch (N = p) has a comparable intensity, which has been 
pointed out only recently (Yakovlev and Levenfish 1995). Rather massive NSs whose neutrino 
luminosity is enhanced by the direct Urea process (Lattimer et al. 1991) 

n — > p + e + u e , p + e — > n + u e . (3) 

belong to the second type. 

At a not too high density p < p , this process is forbidden by momentum conservation. 

Observational data on thermal emission from at least three pulsars - PSR0656 + 14 (Finley et 
al. 1992; Anderson et al, 1993), PSR0833 — 45 (Ogelman and Zimmermann 1989), and Geminga 
(Halpern and Holt 1992; Halpern and Ruderman 1993; Meyer et al. 1994) - indicate that their 
cooling is faster than the standard cooling but slower than that enhanced by the process (|3|). 
In addition, the cooling can be controlled by possible superfluidity of neutrons and/or protons 
in the stellar core. Superfluidity with a moderate critical temperature T c (0 < (T c -T) < T c , 
where T is the temperature of the matter) suppresses the neutrino energy losses but increases 
the heat capacity through the release of latent heat in a phase transition. Strong superfluidity 
(T c T) exponentially suppresses both the neutrino losses and the heat capacity of the NS. It is 
the interplay of these factors that may speed up the standard cooling or slow down the enhanced 
cooling. 

The main parameters of superfluidity that affect the cooling are the critical temperatures 
T cn and T cp of the transition of neutrons and protons, respectively, to a superfluid state. A 
rigorous calculation of T cn and T cp in the core of a NS requires a consistent many-particle quantum 
theory that adequately describes strong particle interactions. Such a theory has not yet been 
constructed, but a large number of theoretical models of superdense matter have been proposed. 
The temperatures T cn and T cp calculated in various models (see Wambach et al. 1991 and Page 
1994 and references therein) vary noticeably with density in the range 10 7 -10 10 K. By comparing 
observations of thermal emission from isolated NSs with calculations of their cooling at various 
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temperatures T cn and T cp , one can place independent constraints on the parameters of neutron 
and proton superfluidity in the NS core. 

A number of studies are devoted to the models of the standard cooling of a NS with a superfluid 
core (see, e.g., Nomoto and Tsuruta 1987 and references therein). Page and Applegate (1992) 
were the first to calculate the cooling enhanced by the direct Urea process in neutron stars with 
a superfluid core. In these studies, the superfluidity of one nucleon species was considered, and 
rather approximate fitting formulas describing the effect of superfluidity on the heat capacity and 
neutrino luminosity were used. A consistent analysis of the effect of n and p superfluidity on 
the heat capacity and neutrino energy losses in the reaction (|3|) was carried out by Levenfish 
and Yakovlev (1993). Gnedin and Yakovlev (1993) used these results and, following Page and 
Applegate (1992), simulated numerically the enhanced cooling of a NS in the presence of neutron 
or proton superfluidity. The enhanced cooling with both n and p superfluidities at once was first 
considered by Van Riper and Lattimer (1993), who based their analysis on very approximate 
factors of suppression of the neutrino luminosity by superfluidity. A consistent calculation of the 
above factors was performed by Levenfish and Yakovlev (1994a,b). Gnedin et al. (1994) used these 
results to model the enhanced cooling of the pulsar PSR0656 + 14 in the presence of a combined 
nucleon superfluidity in its core. Comparison of these calculations with observational data on the 
pulsar's thermal emission allowed the authors to obtain constraints on the critical temperatures 
T cn and T cp in the pulsar's core (assuming that the neutrino energy losses are enhanced by the 
direct Urea process). 

In this paper, we continue to study the cooling of neutron stars with a simultaneous superflu- 
idity of nucleons. Apart from the enhanced cooling, we also consider the standard cooling. The 
suppression factors for neutrino energy losses required for this consideration have recently been 
obtained by Yakovlev and Levenfish (1995). 

Our models are applied to the pulsar 2CG 195 + 04 (Geminga). The effects of neutron and 
proton superfluidities on the standard or enhanced cooling of Geminga were first investigated by 
Page (1994), who used approximate suppression factors for the heat capacity and neutrino energy 
losses. He showed that observations could be explained for any type of cooling by assuming that 
the stellar core was superfluid. We develop more detailed models of the NS cooling by including 
new data on the neutrino energy losses due to the proton branch of the process (0) (Levenfish and 
Yakovlev 1993, 1994a,b; Yakovlev and Levenfish 1995). 

2 MODELS OF COOLING NEUTRON STARS 

We considered two models of a NS. In both cases, the equation of state of matter (npe) in the 
stellar core was taken from Prakash et al. (1988). The maximum mass of the NS for this equation 
of state is 1.7 M . In the first model, the mass of the star is M = 1.44 M , the radius is 11.35 
km, and the central density is p c = 1.37 x 10 15 g-cm -3 . In the second model, M = 1.3 M , 
1Z = 11.71 km, and p c = 1.12 x 10 15 g-cm~ 3 . The chosen equation of state allows the direct Urea 
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process to operate at densities p > p cr = 1.30 x 10 15 g-cm -3 . Thus, the first model corresponds to 
the enhanced cooling: the powerful direct Urea process is possible in a small central kernel with 
the radius of 2.32 km and mass of 0.035 M Q . In the second model, the critical density p cr is not 
achieved, and the star has the standard neutrino luminosity. 

We assumed that nucleons could be superfluid throughout the entire NS core. The proton 
superfluidity results from Cooper pairing of protons in the 1 So state, while the neutron superfluidity 
is due to the neutron pairing in the 3 P2 state with zero projection of the pair's moment onto the 
quantization axis (see, e.g., Wambach et al. 1991). The critical temperatures T cn and T cp can 
depend on density of matter and are not definitely known (Section 1). For certainty, we assumed 
T cn and T cp to be constant throughout the stellar core. 

For our cooling calculations, we used the program (Gnedin and Yakovlev 1993; Gnedin et al. 
1994) in which the inner layers of the NS were assumed to be isothermal. This approximation 
is valid for a not-too-young NS (t > 10-10 3 yr) whose internal thermal relaxation has already 
occurred. Following Glen and Sutherland (1980), the isothermal region was specified by the 
condition p > p b = 10 10 g-cm -3 . 

The photon luminosity of the NS depends on the effective temperature T e of the stellar surface. 
The relation between T e and the temperature Tj near the boundary of the isothermal region 
(p = p h ) is determined by the thermal insulation of the outer envelope (p < pb). The magnetic 
field in the envelope causes the heat conduction to be anisotropic, resulting in an anisotropic 
distribution of T c over the NS surface (see, e.g., Yakovlev and Kaminker 1994). Even a relatively 
weak field B > 10 10 G magnetizes degenerate electrons and reduces the heat conduction across 
B. This enhances the heat insulation of the NS in the vicinity of the magnetic equator. A 
stronger field, B > 10 12 G, is quantizing both for a degenerate electron gas in deep layers of the 
envelope and for nondegenerate sub-photospheric layers. The effects of quantization of electron 
orbits generally increase the longitudinal thermal conductivity of degenerate electrons, reduce the 
pressure of degenerate electrons (in the limit of strong quantization), and enhance the radiative 
thermal conductivity of nondegenerate layers. As a result, the quantizing field weakens the heat 
insulation near the magnetic poles. The NS cooling has been commonly calculated under the 
assumption that the magnetic field is radial in every point of the surface (see, e.g., Van Riper 
1991 and references therein). In this case, only the reduction of the thermal insulation by the 
quantizing magnetic field was taken into account. Recently, Shibanov and Yakovlev (1996) have 
considered the cooling of a NS with a dipole field by taking into account both the reduction of 
the thermal insulation near the magnetic poles and the enhancement of the insulation near the 
equator. Both effects almost compensate each other for a field strength at the magnetic pole 
equal to B ps 3 x 10 12 G, and the cooling proceeds in the same way as for B = 0, although 
the magnetic field gives rise to a strong anisotropy in T e . Our calculations are applied to the 
Geminga pulsar (Section 1) whose magnetic field is close to the above value (see, e.g., Bertsch et 
al. 1992). Therefore, we neglect the effect of the magnetic field and use the relation between 
T c and derived by Van Riper (1991) for B = 0. In this approximation, by T c we mean the 
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average effective surface temperature which determines the total photon luminosity of the star 
L = AirTZ^aTe 4 (without allowing for the gravitational redshift). Note that in the calculations of 
Gnedin et al. (1994), like in other calculations in which the magnetic field has been considered to 
be radial, the influence of the field on the cooling is highly overestimated. 

The neutrino luminosity of the NS with M = 1.44 M & is mainly determined by the direct Urea 
process (|3|). The rate of the corresponding neutrino energy losses in normal (nonsuperfluid) matter 
was calculated by Lattimer et al. (1991). Levenfish and Yakovlev (1993) calculated the factors R 
of suppression of the reaction (|3]) by the neutron superfluidity for normal protons (R = RfP) and 
by the proton superfluidity for normal neutrons (R = R p d ^). They also calculated (Levenfish and 
Yakovlev, 1994b) the suppression of the reaction (|J) by both superfluidities (R = R^). It turned 
out that the reaction is mainly suppressed by the strongest superfluidity, R^ ~ mm[R^\ Rp d ']. 
Note that Van Riper and Lattimer (1993) used the approximate suppression factors R$ ~ R^ ■ 
R p d \ which are especially inaccurate if R^> ~ R p A > <C 1. In addition to the powerful direct Urea 
process, we also took into account the standard neutrino reactions (]!]) and (0). 

The neutrino luminosity of the M = 1.3 M Q star is determined by both the modified Urea 
processes (P and the neutrino emission in nucleon collisions (0) in the entire stellar core. Friman 
and Maxwell (1979) considered the neutron branch of the modified Urea process (the process ([[]) 
for N = n) and the bremsstrahlung emission of nucleons (|2[) in nonsuperfluid matter. Yakovlev 
and Levenfish (1995) studied the proton branch of the modified Urea process (N = p) in non- 
superfluid matter and showed that this branch was almost as efficient as the neutron branch. 
Note that the proton branch was commonly considered to be weak and excluded from calcula- 
tions of the NS cooling. Yakovlev and Levenfish (1995) also calculated the suppression factors 
of the proton (R = _R( m P)) and neutron (R = _R( mn )) branches of the reaction ([[]) by the neu- 
tron superfluidity (R^^ and R^^) for normal protons and by the proton superfluidity (R p mn ^ 
and R^^) for normal neutrons. The suppression factors for the reactions with the same number 
of superfluid particles turned out to be rather close if they were compared as functions of the 
corresponding dimensionless energy gap v = A^/T. Here, A p is the energy gap for the pro- 
ton (singlet) superfluidity, and A n is the minimum energy gap at the Fermi surface of neutrons 
for the anisotropic (triplet) neutron superfluidity. Thus, for example, the factors R^(v p ) and 
R p mn ^ (v p ) are close to each other and to the factors R^(v n ) and R^ p \v n ). Taking this fact into 
account, we propose approximate factors of suppression of the neutron (R = R^ 1 ^) and proton 
(R = R^,^) branches of the modified Urea process by a combined nucleon superfluidity. One can 
expect these factors as functions of the corresponding v (corrected for the number of superfluid 
particles involved) to be not too different from the factor R^J for the direct Urea process. There- 
fore, we use the approximate expressions R^ p \v p ,v n ) ~ R^(2v p , v n ) R^^Vn)/ R$(v n ) and 
R^ n \v p ,v n ) w R$(v p ,2v n ) R p mn \v p )/ RW(v p ). If the protons are normal (v p = 0), the expres- 
sion for R^ p \v p , v n ) becomes exact, while for normal neutrons (v n = 0) the factor R^ n \v p , v n ) 
becomes exact. 

The neutrino energy loss rates in the nucleon-scattering processes (|j) in nonsuperfluid matter 
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are taken from Friman and Maxwell (1979). The loss rate in the pp-scattering is calculated from 
equation (52) of this paper for the nn-scattering with the neutron parameters replaced by the 
proton parameters (see Yakovlev and Levenfish 1995, for details). The factors of suppression of 
the pp- (R = R p pv ^) and np-scattering (R = R^) by the proton superfluidity were calculated 
by Yakovlev and Levenfish (1995). The similarity relations mentioned above allow us to propose 
the approximate suppression factors R^ n ^ for the nn-scattering by the neutron superfluidity and 
the suppression factors R^t^ for the np-scattering by the combined nucleon superfluidity: R^ n ^ ~ 
(v n ) and R^ « R$ (v p , v n ) R^ (v p ) / R^ (v p ) . 

Apart from the neutrino luminosity of the NS core, we take into account the neutrino luminosity 
of the crust resulting from the bremsstrahlung neutrino emission of electrons scattered by atomic 
nuclei. In this case, the approximate formula proposed by Maxwell (1979) is used. 

The heat capacity of the NS is assumed to be the sum of the n, p, and e heat capacities in 
the stellar core; we neglect the heat capacity of the crust because the crust mass is low in the 
chosen NS models. The effects of superfluidity on the n and p heat capacities were investigated by 
Levenfish and Yakovlev (1993). When the superfluidity appears, the heat capacity experiences a 
jump due to the appearance of the latent heat. The proton (singlet superfluidity) and the neutron 
(triplet superfluidity) jump abruptly by factors of 2.43 and 2.19, respectively. At T T c , the 
nucleon heat capacity is exponentially suppressed. In our calculations of the neutrino luminosity 
and heat capacity, we assume the effective masses of neutrons and protons in the NS core to be 
0.7 of the masses of bare particles. 

3 RESULTS 

We calculated about 1800 cooling curves for two NS models (Section 2): with the masses M = 
1.44 M Q (enhanced cooling) and M = 1.3 M (standard cooling). Each curve gives the temporal 
dependence of the effective surface temperature T s = T c ^Jl — TZ g /TZ of the NS as detected by a 
distant observer, !Z g being the gravitational radius of the NS. The curves were calculated in a 
wide range of critical temperatures of neutrons (T cn ) and protons (T cp ) in the NS core (from 10 7 
to 10 10 K). 

Typical cooling curves are plotted in Fig. 1. Curve 2 corresponds to the standard neutrino 
losses (M = 1.3 Mq) in a NS with the normal (nonsuperfluid) core. The change in the shape of the 
curve at t = t u 2.5 x 10 5 yr reflects the change of the cooling regime. For t < t v , the NS is at the 
stage of neutrino cooling: the neutrino luminosity of the inner stellar layers is considerably higher 
than the photon luminosity of its surface. For t > t u , the neutrino luminosity falls well below 
the photon one, and the photon stage sets in. Owing to the low neutrino luminosity, the cooling 
proceeds rather slowly. Curve 1 corresponds to the enhanced neutrino losses (M = 1.44 M Q ) in 
the NS with a normal core. In this case, the NS cools much faster, with the neutrino-cooling stage 
being slightly delayed. The duration t„ of the neutrino stage is mainly determined by the stellar 
heat capacity and by the dependence of the neutrino luminosity on the NS core temperature. For 
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normal nucleons, t„ ~ (2-4) x 10 5 yr. A strong proton superfluidity with normal neutrons reduces 
the heat capacity and time t v by about (20-30)%. For a strong neutron superfluidity and normal 
protons, the heat capacity and t v decrease by a factor of 2 to 3. In the case of strong n and p 
superfluidities, the heat capacity drops approximately by a factor of 20, and t u ~ 10 4 yr. 




lgt(yr) 

Figure 1: Cooling curves for a neutron star with enhanced (curves 1 and 3) and standard (curves 
2 and 4) neutrino luminosities. Curves 1 and 2 correspond to the absence of superfluidity; curve 
4 to T cn = 10 9 - 5 K, T cp = 10 8 - 5 K; and curve 3 to T cn = 10 7 3 K, T cp = 10 8 5 K. The rectangle 
encloses the region of the selected effective surface temperatures T e and Geminga's ages t. 

The results of our calculations are of general character. We will analyze them using the 
Geminga pulsar as an example. Recent observations of Geminga from the ROSAT X-ray Orbital 
Observatory provide a wealth of information on the X-ray emission of this pulsar (Halpern and 
Holt 1992; Halpern and Ruderman 1993; Meyer et al. 1994). The emission consists of soft and 
hard components. The soft component is interpreted as the thermal emission from the Geminga 
surface, while the hard component is commonly attributed to the emission of small hot spots on 



7 



the stellar surface. Theoretical interpretation of the observational data is ambiguous. Thus, for 
example, Halpern and Ruderman (1993) described the soft and hard components by blackbody 
spectra with different temperatures and estimated the effective surface temperature to be T e = 
(5.2 ± 1.0) x 10 5 K. The same authors fitted the soft and hard components with a blackbody 
spectrum and a power-law spectrum, respectively. In this model, T e = (4. 5 ±1.4) x 10 5 K. Meyer et 
al. (1994) interpreted the emission of the hard component using a set of NS hydrogen-atmosphere 
models with a magnetic field. The best-fit models are consistent with a lower temperature, T e = 
(2-3) x 10 5 K. Recently, Page et al. (1995) have interpreted the spectroscopic observations of 
Geminga using the model of a cool hydrogen nonmagnetic atmosphere (T e ~ 2.2 x 10 5 K) without 
a magnetic field in the presence of two small hot spots (4 x 10 5 K) with a radial magnetic field 
of opposite polarity. It should be stressed that the modeling of NS atmospheres is far from being 
complete. It is necessary to solve a complicated problem of radiative transfer in a magnetized, 
sufficiently cool and dense (nonideal) plasma of the NS atmosphere with a complex chemical 
composition (hydrogen, helium, and heavy elements). It is quite possible that new models and 
observational data will yield slightly different values of T e . We will not favor any of the available 
interpretations, but consider the temperature T c as an unknown quantity. We will take T e = (2- 
6) x 10 5 K as the most plausible temperature T s (with allowance for the gravitational redshift). 
For the M = 1.3 M and M = 1.44 M models, the T s /T e ratios are equal to 0.820 and 0.791, 
respectively. 

The age of Geminga is estimated (Bertsch et al. 1992; Page 1994) from the observed rotation 
period P 0.237 s and its time derivative P = 1.10 x 10~ 14 s/s using the well-known formula 
t = (P/P)/(n — 1), where n is the braking index. For the magnetic-dipole losses (n = 3), 
t — ti — 3.4 x 10 5 yr. For n = 4, we have t = t 2 = 2.3 x 10 5 yr, which may be considered as the 
lower limit of the Geming (Page 1994). 

We will assume the ranges of surface temperatures T s = (2-6) x 10 5 K and ages t 2 < t < t\ 
(Fig. 1) to form an "error box" which the real cooling curves of the Geminga pulsar must intersect. 
One can see that in the absence of superfluidity, the cooling proceeds more slowly than required 
for the standard neutrino luminosity and more rapidly than required for the enhanced luminosity. 
A detailed analysis made by Page (1994) has shown that both the standard and enhanced neutrino 
luminosities fit the observations if the effect of superfluidity is taken into account. This assertion 
is illustrated by curves 3 and 4 in Fig. 1. Curve 3 corresponds to the standard neutrino energy 
losses (M = 1.3 Mq) at T cn = 10 95 K and T cp = 10 8 - 5 K. Such a superfluidity decreases strongly 
the heat capacity of the star and accelerates the cooling at the photon stage (when the superfluid 
suppression of the neutrino luminosity is insignificant). Curve 4 corresponds to the enhanced 
neutrino luminosity (M = 1.44 Mq) at T cn = 10 7,3 K and T cp = 10 8 ' 5 K. The superfluidity chosen 
suppresses the neutrino luminosity more strongly than the heat capacity, and delays the cooling 
at the neutrino stage and at the transition stage from the neutrino to photon cooling. 

With the inclusion of new data on neutrino energy losses and the NS heat capacity (Sections 1 
and 2), our calculations extend the study of Page (1994) to a wider range of T cn and T cp . In view 
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Figure 2: T cn and T cp lines corresponding to certain surface temperatures T s of a neutron star 
with the enhanced neutrino luminosity (M = 1.44 M ) and age of 3.4 x 10 5 yr. The numbers near 
the curves are the values of logT s . 



of the great variety of the results, their representation by cooling curves is inconvenient. It is 
more convenient (Gnedin et al. 1994) to display the critical superfluidity temperatures T cn and 
T cp which lead to definite effective surface temperatures T s of a NS of given age t (Figs. 2-5). 
Figures 2 and 4 correspond to the age ti, while Figs. 3 and 5 correspond to the age ti- 

Figures 2 and 3 correspond to enhanced neutrino energy losses. Since the figures are similar, 
we will describe Fig. 2 (t — ti) as an example. The lower left corner of this figure corresponds to 
the case where the critical temperatures T cn and T cp are below or of the same order of magnitude as 
the internal temperature of a NS of age t±. The nucleon superfluidity either has not yet appeared 
in or has had no time to strongly affect the cooling. As a result, the NS surface temperature drops 
to ~ 10 5 K (Fig. 1) through intense neutrino energy losses in the direct Urea process, with the 
star being still at the neutrino cooling stage (t u >t\). 

In all the remaining parts of Fig. 2, the effect of superfluidity on the cooling turns out to 
be crucial. Calculations and a simplified analytic analysis of the cooling equations show that 
t u ~ t\ throughout the figure (with the exception of the upper right corner of Fig. 2). In that 
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case, the cooling is governed by both the neutrino luminosity and the heat capacity of the NS. 
The neutron superfluidity affects the direct Urea process and the proton heat capacity in about 
the same manner as the proton superfluidity affects the direct Urea process and the proton heat 
capacity (Levenfish and Yakovlev 1994a,b). In addition, the heat capacities of normal neutrons 
and protons are comparable but considerably larger than the electron heat capacity. As a result, 
the curves in Fig. 2 are almost symmetric with respect to the replacement of the coordinate axes, 
logT cra <-> logT cp , and it is sufficient, for instance, to describe the lower right half of Fig. 2. 

If a low temperature of T cp ~ 10 7 K is fixed and T cn increases from ~ 10 8 to ~ 10 10 K, the 
proton superfluidity does not occur by the time t±, and the cooling is governed by the neutron 
superfluidity alone. Therefore, the T s = const lines in the lower right corner of Fig. 2 are vertical. 
As T cn increases, the neutrino luminosity drops by several orders of magnitude. The heat capacity 
also changes but much weaker, because the strong n superfluidity reduces the NS heat capacity 
only by a factor of approximately 3. As a result, variations in the neutrino luminosity are more 
pronounced than variations in the heat capacity. The cooling of the star slows down, while T s 
rises. If, on the other hand, we fix T cn > 10 8 K and increase T cp , then the p superfluidity will arise 
in a star of age t\ at some value of T cp (and the T s = const line will no longer be straight). As 
T cp increases further, the superfluidity suppresses the neutrino luminosity and influences the heat 
capacity. However, the additional suppression of the direct Urea process by a moderate proton 
superfluidity in the presence of a strong neutron superfluidity is weaker than the change in the 
proton heat capacity (Levenfish and Yakovlev 1994a,b). Then the cooling is mainly controlled by 
the proton heat capacity. When the superfluidity arises, the heat capacity increases in a jump-like 
manner (the latent heat is liberated), giving rise to an increase in T s : the T s = const lines are 
shifted leftward. However, the increase in the heat capacity is superseded by its sharp suppression 
with rising T cp . The NS becomes cooler, and the isolines are shifted rightward. 

Finally, the critical temperatures T cn and T cp in the upper right corner of Fig. 2 are so large 
that the neutron and proton superfluidities suppress the direct Urea process and the nucleon heat 
capacity almost completely. For a very low (electron) heat capacity, the neutrino-cooling stage is 
short (t u ~ 10 4 yr, see above). As a result, the NS had long cooled down by the time ti through 
the photon emission, and its temperature is very low, T s ~ 10 5 K. 

Figures 2 and 3, corresponding to the enhanced cooling of Geminga, are similar to Figs. 2, 4, 
and 6 from Gnedin et at. (1994) composed for the younger pulsar PSR0656 + 14 (t ks 10 5 yr). A 
comparison of these figures shows that higher critical temperatures of one of the nucleon compo- 
nents under the condition T cp <C T cn or T cn <C T cp are required for maintaining the high effective 
temperature T s over a longer period of time. 

T s isolines for the NS model with M = 1.3 M Q are plotted in Figs. 4 and 5. Figures 4 and 
5 correspond to the ages t\ and £2, respectively. It can be seen that the behavior of the isolines 
for the standard neutrino energy losses differs qualitatively from that for the enhanced neutrino 
energy losses (Figs. 2 and 3). The principal difference lies in the fact that the relatively weak 
standard neutrino luminosity of a NS of age t\ or t 2 turns out to be lower than the photon 
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luminosity (t„ < t 2 ). The NS of the above age is at the photon-cooling stage virtually for all 
values of T cn and T cp shown in Figs. 4 and 5, and the neutrino luminosity of the stellar core plays 
no significant role. The nucleon superfluidity affects the cooling mainly through the heat capacity. 
The exception is only some small parts of Figs. 4 and 5. 

Figures 4 and 5 are qualitatively similar. We will describe Fig. 4 as an example. At T cn < 
10 7 - 9 K and T cp < 10 7 9 K (lower left corner of Fig. 4), the nucleon superfluidity in the core of a 
NS of age t-y has not yet appeared, and the surface temperature is independent of T cn and T cp . At 
T cp ^ 10 7 ' 5 K (lower part of Fig. 4), the protons are normal. If, at the same time, we raise T cn , the 
neutron superfluidity will arise at T cn > 10 79 K by the time t\. At T cn < 10 8 K, this superfluidity 
is moderate. The latent heat is released, and T s slightly rises. As T cn increases further, the n 
superfluidity becomes strong by the time t\\ it sharply suppresses the neutron heat capacity, and 
decreases T s . As T cn rises, the n superfluidity arises increasingly earlier. At T cn ~ 10 9 K, it 
appears at the early neutrino-cooling stage and is able to strongly suppress the neutrino losses 
and delay the cooling. Though weakly, this delay manifests itself by the time t\. as T cn becomes 
higher than 10 9 K, the temperature T s slightly rises. If, however, we fix T cn ~ 10 9 K and raise 
T cp , the p superfluidity arises at T cp > 10 7 6 K. The latent heat is released, and T s grows. As a 
result, the T s = const lines in Fig. 4 form loops (e.g., logT s = 5.7). 
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At T cn < 10 7 ' 9 K (left part of Fig. 4), the neutrons are normal. If T cp is raised above > 10 7 ' 9 K, 
the proton superfluidity appears by the time t\. At T cp ~ 10 s K, the superfluidity arises just before 
t\ and is moderate. The heat capacity increases due to the latent-heat release, and T s slightly rises. 
As T cp grows further, the p superfluidity arises at increasingly earlier stages and becomes strong 
by the time t±. It suppresses the proton heat capacity, which amounts to < 30% of the total 
heat capacity of the NS for normal neutrons. On the other hand, this superfluidity exponentially 
suppresses the neutrino luminosity at the preceding neutrino-cooling stage. The slowdown of the 
cooling due to the suppression of the neutrino luminosity is more pronounced than the acceleration 
caused by an insignificant suppression of the heat capacity, and T s rises with increasing T cp . 

At T cp > 10 9 K (upper part of Fig. 4), the p heat capacity is suppressed completely. If T cn is 
raised, the n superfluidity appears by the time t±, which first enhances and then suppresses the 
n heat capacity (see above). This causes T s to initially rise and then decrease with a peak at 
T cn « 10 8 - 3 K. If T cn « 10 8 - 3 K is fixed and T cp decreases, T s slightly drops. This drop of T s is 
caused by the enhancement of the stellar neutrino luminosity at the preceding neutrino-cooling 
stage. As a result, the logT, = 5.9 curve forms, for example, a loop in the upper part of Fig. 4. 
The logT s = 5.77 loop in the right part of Fig. 4 is explained in a similar way (with replacement 
n <-> p). 
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Finally, at T cn > 10 8 6 K and T cp > 10 8 ' 3 K (upper right part of Fig. 4), the n and p su- 
perfluidities become strong by the time t\ and heavily suppress the nucleon heat capacity. The 
neutrino-cooling stage is considerably shorter {t v <C t±, see above), and the T s = const curves are 
entirely determined by a decrease in the heat capacity caused by superfluidity. The higher T c , the 
lower the heat capacity, and the cooler the star. At T cn > 10 92 K, the n superfluidity suppresses 
the neutron heat capacity almost completely, and the cooling is governed only by a change in T cp 
(proton heat capacity) at T cp < 10 9 K. If, however, T cp > 10 9 K, the p superfluidity suppresses the 
proton heat capacity almost completely, and the cooling is governed by a change in T cn (neutron 
heat capacity) at T cn < 10 92 K. If the combined superfluidity of neutrons and protons is very 
strong (T cn > 10 9 ' 2 K, T cp > 10 91 K, see the uppermost right corner of Fig. 4), the nucleon heat 
capacity of the NS core is fully suppressed, and the cooling no longer depends on T cn and T cp . 
In this regime, the cooling is governed by the electron heat capacity, which is independent of the 
nucleon superfluidity. 

It should be emphasized that the T s = const curves at T cn > 10 8 - 5 K and T cp > 10 8 3 K (upper 
right part of Fig. 4) are determined by an abrupt suppression of the heat capacity by strong 
superfluidity and are only slightly sensitive to the specific parameters of the NS. In the remaining 
region of T cn and T cp , variations in T s are small, but the T s = const curves (for instance, loops) 
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may formally strongly deform with changing NS parameters (cf., the logT s = 5.7 and logT, = 5.8 
isolines in Figs. 4 and 5). This deformation is not accompanied by large variations in T s and does 
not affect the main conclusions concerning the effect of superfluidity on the cooling. 

Note that Page (1994) considered a wider age range for Geminga, t 2 < t < t 3 , where t 3 = 
6.8 x 10 5 yr corresponds to the braking index n = 2 (see above). We give no results for t = t 3 . In 
the case of the enhanced neutrino losses, the corresponding figures are similar to Figs. 2 and 3. 
In the case of the standard energy losses, at t = t 3 , Geminga is at the advanced photon-cooling 
stage, when the surface temperature is sufficiently low (T s « 5.7 x 10 5 K) even for normal nucleons 
(Fig. 1). 

4 CONCLUSIONS 

Our calculations support the conclusion of Page (1994): observational data on the thermal emission 
of Geminga can be explained by both the standard or enhanced neutrino energy losses if neutrons 
and protons in the Geminga's core are assumed to be superfluid. Figures 2 and 3 give the critical 
temperatures of the neutron and proton superfluidities, T cn and T cp , which lead to specified effective 
surface temperatures T s of Geminga if its age is ti = 3.4 x 10 5 yr or t 2 = 2.3 x 10 5 yr (Section 3) 
and the neutrino luminosity is enhanced by the direct Urea process (|j). Figures 4 and 5 are similar 
but correspond to the assumption that the neutrino luminosity is standard and results from the 
reactions ([!]) and (H). Observational data on the X-ray emission of Geminga (Section 3) indicate 
that its effective surface temperature T s ranges from 2 x 10 5 K to 6 x 10 5 K. 

Both assumptions about the neutrino energy losses (standard and enhanced) can be reconciled 
with observational data. In both cases, the required values of T cp and T cn are in qualitative 
agreement with various model calculations of the critical temperatures of nucleon superfluidity in 
NS cores (see Page and Applegate 1992, for references). However, the properties of superfluidity 
in the Geminga's core for the standard and enhanced neutrino luminosities must greatly differ. 
In the case of the enhanced neutrino luminosity, the superfluidity is needed to slow down the 
NS cooling at the neutrino-cooling stage (or at the transition stage from the neutrino to photon 
luminosity) through a partial suppression of the neutrino energy losses in the direct Urea process. 
For example, the values of T cn and T cp that give the Geminga's surface temperature T s = 10 5,3 - 
10 5 ' 5 K lie (see Figs. 2 and 3) in two relatively wide decoupled regions in the (T cn -T cp ) plane. 
For the standard neutrino processes, the superfluidity must slow down the cooling at the photon- 
cooling stage through the suppression of the stellar heat capacity. The region of the necessary 
values of T cn and T cp (Figs. 4 and 5) turns out to be much narrower and coupled. For the same 
T s , we obtain logT cn > 8.7 and logT cp > 8.5. In this case, the constraints on possible values of 
critical temperatures are much more stringent. On the whole, the model of the standard neutrino 
losses requires higher T cn and T cp than the model of the enhanced losses. 

The results of our calculations are tentative. More detailed models for the cooling of NSs with 
superfluid cores are needed. First, it would be desirable to allow for the nonuniform distribution 
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of the critical temperatures T cn and T cp (for their density dependence) throughout the stellar 
core. For this to be done, joint theoretical calculations of T cn and T cp are necessary for various 
models of superdense matter (which are now not available). These data could be used to directly 
determine the best-fit models without assuming that T cn and T cp are free parameters. It would be 
desirable to consider the cases where other particles, above all, muons and hyperons, are present 
in the NS cores along with neutrons, protons and electrons. Note that the neutrino-loss rates in 
various types of reaction are also sensitive to superdense-matter models and have not yet been 
determined with confidence. At the late cooling stage, additional heating sources may also be 
important, for example, those related to dissipation of the rotational energy in the core of a star 
during its slowdown (see, e.g., Cheng et al. 1992, Umeda et al. 1994, and references therein). 
One may hope that new observational data on the thermal emission of NSs together with new 
theoretical calculations will make it possible to obtain more stringent constraints on the rate of 
neutrino energy losses and the critical temperatures of nucleon superfluidity in superdense matter. 
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